Thermodynamics of Coarse Grained Models of Super-Cooled Liquids 
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In recent papers, we have argued that kinetically constrained coarse grained models can be applied 
to understand dynamic properties of glass forming materials, and we have used this approach in 
various applications that appear to validate this view. In one such paper [J. P. Garrahan and D. 
Chandler, Proc. Nat. Acad. Sci. USA 100, 9710 (2003)], among other things we argued that 
this approach also explains why the heat capacity discontinuity at the glass transition is generally 
larger for fragile materials than for strong materials. In the preceding article, Biroli, Bouchaud and 
Tarjus (BB&T) have objected to our explanation on this point, arguing that the class of models we 
apply is inconsistent with both the absolute size and temperature dependence of the experimental 
specific heat. Their argument, however, neglects parameters associated with the coarse graining. 
Accounting for these parameters, we show here that our treatment of dynamics is not inconsistent 
with heat capacity discontinuities. 



I. INTRODUCTION 

Many workers argue that a thermodynamic anomaly 
underlies the onset of glassy dynamics |]|. Evidence for 
this view is the rough, though not quantitative |2| , corre- 
lation between dynamic fragility and excess heat capac- 
ity discontinuity |3| ■ This thermodynamic view contrasts 
with the picture we have advocated 0, attributing 
glassy behavior to dynamic heterogeneity |4j , with grow- 
ing length scales appearing in space-time, but not space 
alone |5| . Indeed, direct observations of diffusive motion 
in colloidal glasses reveal excitations that are local and 
sparse 7j. These findings seem consistent with excita- 
tions being local and uncorrelated, as assumed in the two- 
state model for the low temperature behavior of struc- 
tural glass heat capacities @- Spatial correlations under 
such conditions can arise through constraints on particle 
motions that are relieved only when adjacent regions have 
exhibited some degree of mobility 0)ll3- Trajectories are 
then correlated throughout space and time, with varying 
degrees of hierarchical structure determining the extent 
to which the system is fragile or strong From this 
perspective, a non-thermodynamic explanation emerges 
for correlation between heat capacity and fragility |12| : 
The concentration of excitations required to support frag- 
ile hierarchical dynamics is higher than that required for 
strong non-hierarchical dynamics. The juxtaposition of 
heat capacity discontinuities is therefore understood as a 
consequence of different excitation concentrations. 

Biroli, Bouchaud and Tarjus (BB&T) |l3| have taken 
issue with this explanation. They show that the simplest 
treatment of defect models cannot simultaneously fit the 
size and temperature dependence of the experimental 
heat capacity while simultaneously fitting dynamic prop- 
erties. In general and in more current context, however, 
the terminology "defect models" refers to a broad class 
of kinetically constrained models [l4| . As a result of this 
generality, their question "Are defect models consistent 
with the entropy and specific heat of glass-formers?" has 
a non-trivial answer (see also the analysis in (lsHIij. Il7| ) . 



In particular, this class of models has an assortment of 
possible dynamical behaviors, so that relaxation data can 
be fit with many different functions of excitation concen- 
tration. In addition, kinetically constrained models are 
coarse grained caricatures of fluids, so that many degrees 
of freedom remain unspecified. Here, we focus on this 
latter feature, the possible consequences of unspecified 
degrees of freedom. 

The starting point is the assumption that small length- 
scale and small time-scale features of a fluid can be in- 
tegrated out leaving only simple stochastic rules for the 
dynamics of discrete variables on a lattice. These rules 
contain constraints imagined to be the consequence of the 
actual intermolecular interactions, interactions that limit 
the space or metric for molecular motions 0] . The exci- 
tations or defects that survive coarse graining distinguish 
microscopic regions of space-time that exhibit molecular 
mobility from those where molecules are jammed or im- 
mobile. In particular, an excitation or defect is a micro- 
scopic region of space for which particle mobility emerges. 
This characterization is related to coarse graining in time, 
because it takes time to discern whether or not mobility 
is exhibited. Bear in mind, "defect" does not distinguish 
disordered arrangements of atoms from those that are 
ordered because there are generally many disordered yet 
jammed configurations. Similarly, a lack of "excitation" 
does not necessarily imply low energy because some im- 
mobile regions may have the same energy as mobile re- 
gions. Uncertainty concerning this terminology may be 
the ori gin of criticisms of facilitated models leveled by 
BB&T [TUTU and also by Lubchenko and Wolynes [2(f . 

As a simple illustration, consider dividing space into 
cells, with grid spacing larger than the equilibrium cor- 
relation length of the material. At a given time frame, 
the micro states of different cells are then uncorrelated. 
Further, suppose there are two energy levels, and J, 
for the states of a given cell, where the states with en- 
ergy J include those that will exhibit particle mobility 
after a coarse graining time St. There is a fraction, 0, of 
states with energy J that exhibit mobility in this way; 
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the others remain jammed. In this case the average exci- 
tation concentration is <f> times the concentration of cells 
with energy J, and thus the number of molecules con- 
tributing to energy or enthalpy fluctuations is times 
larger than that contributing to mobility fluctuations. In 
other words, when considering the thermodynamic impli- 
cations of a kinetically constrained defect model, a factor 
<f>~ 1 is required to account for energetic states or degrees 
of freedom absent from explicit consideration. BB&T ig- 
nore this factor, which leads them to the difficulties they 
describe. 

In the next section we discuss this point in greater de- 
tail. We show that facilitated or kinetically constrained 
models such as the one we presented in Ref. ^2 a d- 
mit a variety of possibilities for partitioning micro states. 
This demonstration justifies the heat capacity formulas 
we used in Ref. 0] and therefore supports our argument 
about the juxtaposition of heat capacities for strong and 
fragile glass formers. 



While the state \x of a cell will correspond to one of 
the two mobility states, n = 0,1, in general its energy, 
tfj_, will be distributed. (For notational simplicity, when 
clarity is not diminished, we drop the subscripts on rii 
and fii.) Let Go(e) and Gi(e) represent the normalized 
probability densities for the energy levels in states n = 
and n = 1, respectively. The cell partition function then 
reads (we set fcs = 1): 



Z = J2^ /T = z Q + z u 



where 



de G n (e)e 



e/T 



(n = 0,l), 



(1) 



(2) 



and w n gives the total number of levels in each of the two 
mobility states. The excitation concentration is given by: 



II. MODELS 

The dynamics of facilitated models we consider are 
governed by master equations for distribution functions 
of the mobility field on a lattice. Consistent with our ear- 
lier papers, we use the symbol riiit) to denote the value 
of the mobility field at lattice site i during time period t. 
It is a binary field in that rii takes on one of two values, 

(corresponding to a cell that is unexcited, jammed) or 

1 (corresponding to excited, mobile). Each cell is imag- 
ined to contain several molecules, so we expect that even 
after some coarse graining, there are still many (not just 
two) microstates possible for each cell at a given time 
t. We further imagine these many microstates, specified 
with /Xj, evolve according to their own master equation. 
The master equation for the distribution of the mobility 
field rii is a contraction of that for the micro-state field 
//j. In this section, we show how this contraction works. 
We begin by considering the thermodynamics and statis- 
tics of the fields rii and fii, and then consider the master 
equations for the distributions of these fields. 



A. Excitations and microstates 

Imagine we partition the liquid into a lattice with grid 
spacing 5x, and denote the state in cell i by fii. Imagine 
further that the set of available cell levels {fii} can be 
split into two groups: one subset of {fii} is "unexcited" 
with respect to mobility, m = 0, and the complementary 
subset is "excited", m = 1. In practice, rii is most con- 
veniently determined by observing system for a coarse 
graining time St. This splitting, or projection, of many 
micro states into two different mobility states (or several, 
as in the generalization of is the the first central as- 
sumption of facilitated models 0] . 



(n) 



Zl/ZQ 



1 + Zl/Zo' 

The average energy per cell is: 

(e) = (l-c)(e )+c<ei> 

where 

/k ,_Jde e k G n {e)e~^ 



(n = 0,l). 



(3) 



(4) 



(5) 



xn '~ J de G n (e)e-*/ T 
The specific heat per cell is: 

C v = T~ 2 « £1 ) - <6o)) 2 C (l- C ) + C(°)(l-c) + C«c, (6) 
where 

T 2 C^ ee (4) - ( e „) 2 (n = 0,l). (7) 

The quantity ]h(zi/zq) determines the concentration 
c of excitations; (ei) — (eo) gives the average cell energy; 
and T- 2 ((ei) - (e )) 2 + C ( v 1] - C { ° ] sets the specific heat 
at small c. Except in the case where G n (e) are delta 
functions, these scales are in principle all different. 

As a simple example, which will be useful below, con- 
sider the following model with just two energy levels. The 
unexcited states n — occupy both energy levels, which 
we denote g and h, with energies e g = and e/j = J, 
respectively: 



G (e) = (l-a)S(e)+aS(e-J), 



(8) 



where a sets the relative degeneracy between the two 
levels. The excited states n = 1 have a single energy 
level, which we denote x, with energy e x = J: 



Gi(e) = S(e - J). 



(9) 
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In this case, the cell partition function, average energy 
and specific heat become: 



Z = (1 — a)ujQ + e ' ] I T (aujQ + w{) 
(e) = J(l-c)c h + Jc 
C v = {J/Tf{l-c h fc(l-c) 

+ {J/T) 2 c h (1 - Ch) (1 - c), 



(10) 
(11) 

(12) 



where cu denotes the relative occupation of the states h 
with respect to states g, Ch = ae~ J / T /[(l — a) +ae~ J ' T ]. 
The total concentration of levels with energy J is 



■+ (1 - c)c h = <f> 1 c, 



(13) 



where we have defined cj> as the ratio of mobile cells to 
the total number of cells with energy J. The specific heat 
then reads 



C v = {J/Tf 0- 1 c (1 - (T 1 c). 



(14) 



Kinetics 



We now show how an appropriate projection of the 
dynamics of the full states /i, of individual cells into 
their mobility states n, = 0, 1 gives a two-state facilitated 
model such as the FA or East models, and the slightly 
more general class of models introduced in Ref. 12] . The 
approach we take has been used to derive approximate 
real-space renormalization group equations |22j. Here, 
however, the approach requires no approximation. 

The master equation for the dynamics can be written 

m. 



d_ 

Of 



P(t)) = -C\P(t)), 



(15) 



where \P(i)) is the state vector for the probability density 
of the system at time t, and C is the Liouvillian operator 
for the dynamics. In the direct product representation it 
reads 22]: 

Mi)®IM2)®- • -<8>|/xjv> , (16) 

M 

where is the state vector for level fa in site i. 

The second assumption in facilitated models is that 
transitions in a cell i are only possible if neighbouring 
cells j are mobile, n,j = 1 |14l . In a system with facili- 
tated dynamics C reads I22j: 



(17) 



where Ci is the kinetic constraint operator imposing fa- 
cilitation, and £i is the unconstrained dynamic opera- 
tor at site i. For example, in the East facilitated model 
Ci = rii-i, and in the direct product representation: 



(g) m 



(18) 



The off diagonal elements of li are given by the tran- 
sition rates between levels: ili)^^ = —7^^^', while 
the diagonal elements are such that the sum of columns 
in £i vanish. A choice satisfying detailed balance is 
(^i)p'm = ~ (A 4 ' l ec l): where |eq) is the equilibrium dis- 
tribution. Notice that we have assumed that transitions 
between any of the levels in the cell are only possible if 
facilitated by a neighbour. 

In order to project the dynamics of the full set of levels 
into that of only the excitation occupancies \n = 0, 1) 
we need projection, T\, and embedding, T2, operators 

M 

T i=^2 a n,ix W) (^1 > T 2 = ^2 b^.„ \fi) (n\ , (19) 



71, fl 



11 . n 



where the a and 6's are constants. The projection is 
performed independently in each cell, so these are site 
diagonal operators. They must also satisfy TiT 2 = 1. 
The dynamics of the projected state \P(t))' = T\ \P(t)) 
is then given by the projected Liouvillian 



C! = Ti C T 2 , 



(20) 



where the elements of the projection and embedding ma- 
trices ifT^jl are given by: 



(/u|eq) 



l-n 
1 -c 



(21) 



The choice of rates in C, which depend only on the fi- 
nal level, j^fj,' — (// |eq ), ensures that the projected 
dynamics is Markovian 23]. 

We now show explicitly how the projection procedure 
works in the simple the three-level example of the previ- 
ous subsection. In this case = \g) , \h) , |x), and the 
onsite operators can be represented by 3 x 3 matrices: 



(22) 




£ = 



c+(l-c)c h -(l-c)(l-c h ) 

-(1-C)c h C+(1-C)(l-Cfc) 

— c — c 

-(l-c)(l- Ch ) 
-(1 - c)c h 
l-c 



(23) 



Of the three levels only \x) corresponds to the 
excited state, and the projection matrix is: 



Ti 



1 1 
1 



The corresponding embedding matrix is: 



T 9 



1-Cfc 
c h 



(24) 



(25) 
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The projected operators then read: 
n' = T\nTi ■■ 
£' = T X IT 2 = 




1 



c c — 1 
— c 1 — c 



(26) 
(27) 



These are precisely the operators for a two state facili- 
tated model IH. 



III. DATA FITTING 

Subsection III Al tells us that if the mobility field is the 
result of coarse-graining and projecting the many mo- 
bile and immobile micro states, then different energy (or 
enthalpy) scales enter into the concentration of mobile 
cells c, the average energy (or enthalpy), and the spe- 
cific heat. How these scales relate to each other depends 
on the distribution of energy levels. At equilibrium, the 
average excitation concentration, c, is a function of ther- 
modynamic state. Its value should be indicative of the 
concentration of regions with excess energy or enthalpy, 
but to say more requires more assumptions or detail than 
needed to predict dynamical trends. 

In order to illustrate how facilitated models can be 
used to fit thermodynamic data the simple example given 
above of two immobile cell levels g and h and a single mo- 
bile level x will suffice. From Eqs. (| 1 Of) — (| 1 4|) we see that 
while c gives the concentration of cells displaying mo- 
bility, both mobile and immobile cells can contribute to 
energy (or enthalpy) fluctuations. The relevant quantity 
to determine the specific heat is not c but c + (1 — c)ch- 
From Eqs. and (|14ft we have that the specific heat 
per cell is: 



Cecil « f^J/rfcfl-f'c) 

~ r 1 {j m 2 c, 



(28) 



The asymptotic equality of Eq. (|28(1 is, in effect, the ex- 
pression we refer to in the penultimate paragraph of Ref. 
[l2| . The factor _1 accounts for the number of molecules 
contributing to energy or enthalpy fluctuations, this num- 
ber being larger than that for molecules contributing to 
mobility fluctuations. This fact is implicit in the coarse 
graining procedure, as shown above. In Ref. |T^. is 
included in the quantity we called A/". In the few lines 
devoted to this topic in Ref. 0] , the meaning of N was 
obscure and its description misleading. Indeed, where we 
wrote "A/" is the number of molecules that contribute to 
enthalpy fluctuations per mobile cell", we should have 
written "A/" accounts for the number of molecules that 
contribute to enthalpy fluctuations beyond those associ- 
ated with mobile cells." 

To be perfectly explicit, if s is the average number of 
molecules in a cell, then N — l/(j>s. With this expression, 
the relative values of C for several super-cooled liquids 
at T g are consistent with the relative values of their heat 
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FIG. l: Fitting ol tnermoayn 

facilitated models of Ref. fl2|| and compared to the formula 
used by Biroli, Bouchaud and Tarjus (indicated as BB&T) 
|l3| . Top panel: excess specific heat, Ac p , of supercooled 
OTP (left hand scale) and 3BP (right hand scale). Symbols 
indicate experimental data lUI . Full curves are fits using the 
formulas for the heat capacity per molecule discussed in the 

A Cp ): c mo i. = A/y/r) VMi - <T lc ), 



and 



<j> e + 4>i(T-T s )/T g . Here 



text (where C m oi. 
where N = (4>s)' 

c w (j}e~ J/T uj 1 /[(l -a)u ] ~ g exp (- J/T + J/T Icf ), where 
the first (approximate) equality neglects terms of order c 2 , 
and the second introduces the notation used in our Ref. [T^|: 
g is a degeneracy that we imagine measures the possible di- 
rections of mobility, and T rc f bounds from above the tem- 
peratures we consider. The values of J « 26.7 and 16.7 for 
OTP and 3BP were determined in |T^| from kinetic data. On 
physical grounds we expect s « (J/T m )/AS , f usioll ~ 1 — 10, 
and (j> to increase with temperature. The fitting of Ac p gives 
s, 4> g ,4>i = 5, .03, .7 and 4.5, .04, .5 for OTP and 3BP, respec- 
tively. JV ~ 6 at T g as argued in 0. Clearly, it is pos- 
sible to describe the size and modest change with T of the 
specific heat in the vicinity of T g , in contrast to the observa- 
tion of BB&T. Bottom panel: excess entropy, S c , of OTP 
(left hand scale) and 3BP (right hand scale). Symbols cor- 
respond to experimental data 16]. Full curves are fits us- 
ing the formula for the entropy per molecule which follows 
fromEq. J2HJ: S mo l. = s _1 {ln[(l - a)u ] - (^c) In (0 _1 c) - 
(1 - ^c) In (1 - (P^c) + (<^- 1 c) In [(w a + wi)/(l - a)w ]}- 
The shown fits correspond to the same parameters as be- 
fore plus s _1 In [(1 — a)coo], s _1 In [(u/oa + wi)/(l ~ a)^o] = 
2.8, 1.2 and 2.4, 1.8 for OTP and 3BP, respectively. 
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capacity jumps at the glass transition. This consistency 
is independent of the value of s and Af, but the value of 
Af does play a role in the actual value of the heat capacity 
discontinuity. Specifically, 

C mol . « a^Ccdi ~Af{J/T) 2 c, (29) 

is consistent with the value of the discontinuities of the 
heat capacity per molecule at T g when Af w 6. BB&T 
[l3| identify s with A/", and further neglect the factor 1/0, 
and in so doing disagree with this consistency. 

BB&T's second criticism is about the temperature 
variation of C. They consider the relationship 

Cbb&t tx ( J/T) 2 c. (30) 

This formula is consistent with the asymptotic equality 
of Eq. J2HJ, which requires 4>~ 1 c to be small. While it 
is small at the glass transition, <p~ 1 c = scAf may not be 
small at the higher temperatures considered. The right- 
hand side of BB&T's formula is missing the requisite fac- 
tor of (1 — _1 c) present in Eq. It is also missing a 
possible temperature dependence of <f>. This temperature 
dependence would reflect that of transport in a liquid. 
It is typically more pronounced for isobaric temperature 
variation than for isochoric temperature variation, but 
weak and sub-Arrhenius in either case [2J|. Thus, we 



write <p w g + 4>\ (T — T s ) /T g , where the parameters <f> g 
and <pi are temperature independent. 

In Ref. ^3 j we chose to not consider C above T g be- 
cause doing so requires the additional parameters de- 
scribed in the previous paragraph. But since, BB&T 
raise this point, we show in Fig.l that the effects of 
these extra factors on the temperature dependence of C 
is significant. This figure can be contrasted with that in 
BB&T's ^3|. While the asymptotic formula predicts an 
exponential variation with temperature, the full formula 
predicts modest variation with temperature that is not 
inconsistent with experiment. We also show fits to the 
corresponding excess entropy, S c . Better agreement be- 
tween theory and experiment for these properties would 
require modification of the models and perhaps further 
parameterization. 



Acknowledgments 

In the US, this work was supported initially by the 
NSF, and more recently by DOE grant no. DE-FE-FG03- 
87ER13793. In the UK, it was supported by EPSRC 
grants no. GR/R83712/01 and GR/S54074/01, and Uni- 
versity of Nottingham grant no. FEF 3024. 



[1] See, for example, D. Kivelson, S.A. Kivelson, X.-L. Zhao, 

Z. Nussinov, and G. Tarjus, Physica A 219, 27 (1995); M. [10 

Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999); 

X. Xia and PG. Wolynes, Proc. Natl. Acad. Sci. USA [11 

97, 2990 (2000); S. Franz and G. Parisi, J. Phys. C 12, 

6335 (2000). [12 

[2] D. Huang and G.B. McKenna, J. Chem. Phys. 114, 5621 

(2001) ; H. Tanaka, Phys. Rev. Lett. 90, 055701 (2003). [13 
[3] See e.g. C.A. Angell, Science 267, 1924 (1995). 

[4] For reviews see: H. Sillescu, J. Non-Cryst. Solids 243, 

81 (1999); M.D. Ediger, Annu. Rev. Phys. Chem. 51, [14 
99 (2000); S.C. Glotzer, J. Non-Cryst. Solids, 274, 342 
(2000); R. Richert, J. Phys. Condens. Matter 14, R703 [15 

(2002) . [16 
[5] J. P. Garrahan and D. Chandler, Phys. Rev. Lett. 89, 

035704 (2002). [17 
[6] L. Berthier and J.P Garrahan, Phys. Rev. E 68, 041201 [18 

(2003) ; S. Whitelam, L. Berthier and J.P. Garrahan, 
Phys. Rev. Lett. 92, 185705 (2004); Y.J. Jung, J.P. Gar- [19 
rahan and D. Chandler, Phys. Rev. E 69, 061205 (2004); 

L. Berthier, D. Chandler and J.P. Garrahan, Europhys. [20 

Lett. 69, 230 (2005). 
[7] E. Weeks, J.C. Crocker, A.C. Levitt, A. Schofield, and [21 

D.A. Weitz, Science 287, 627 (2000). [22 
[8] P.W. Anderson, B.I. Halperin and C. M. Varma, Philos. 

Mag. 25, 1 (1972); W.A. Phillips, J. Low Temp. Phys. 7, [23 

351 (1972). 

[9] S.H. Glarum, J. Chem. Phys. 33, 639 (1960); M.C. [24 
Phillips, A.J. Barlow and J. Lamb, Proc. R. Soc. Lond. 



A 329, 193 (1972). 

G.H. Fredrickson and H.C. Andersen, Phys. Rev. Lett. 
53, 1244 (1984). 

R.G. Palmer, D.L. Stein, E. Abrahams and P.W. Ander- 
son, Phys. Rev. Lett. 53, 958 (1984). 
J.P. Garrahan and D. Chandler, Proc. Natl. Acad. Sci. 
USA 100, 9710 (2003). 

G. Biroli, J.-P . Bouchaud and G. Tarjus, 
cond-mat/0412024 and J. Chem. Phys. preceding 
article. 

For a review see F. Ritort and P. Sollich, Adv. Phys. 52, 
219 (2003). 

C.A. Angell and K.J. Rao, J. Chem. Phys. 57, 470 (1972). 

C. T. Moynihan and C.A. Angell, J. Non-Cryst. Solids 
274, 131 (2000). 

D. V. Matyushov and C.A. Angell, |cond-mat /0503082 

S. Whitelam and J. P. Garrahan, J. Phys. Chem. B 108, 
6611 (2004). 

G. Biroli and J.-P. Bouchaud, J. Chem. Phys. 121 , 7347 
(2004). 

V. Lubchenko and PG Wolynes, J. Chem. Phys. 121 , 
2852 (2004). 

E. D. Siggia, Phys. Rev. B, 16, 2319 (1977). 

See, for example, S. Whitelam and J.P. Garrahan, Phys. 
Rev. E 70, 046129 (2004). 

M. Sharpe, "General Theory of Markov Processes" (Aca- 
demic Press, Boston, 1988). 
J. Jonas, Science 216, 1179 (1982). 



